Dynamic structure factor of Fermi superfluid in the BEC-BCS crossover 
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We consider cigar shaped Fermi superfluid in the BEC-BCS crossover. Using polytropic form of 
equation of state, we derive low energy multibranch bosonic excitations and the corresponding 
density fluctuations in three different regimes along the crossover, namely weak-coupling BCS, 
unitarity and molecular BEC regimes. Bragg spectroscopy can be used to probe the multibranch 
nature of the low energy bosonic excitations by measuring dynamic structure factor. Therefore, we 
calculate dynamic structure factor in those three different regimes. In the Bragg spectroscopy, an 
actual observable is momentum imparted to the superfluid due to the Bragg potential. We also 
present results of the momentum imparted to the superfluid due to the Bragg pulses. 
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I. INTRODUCTION 

The crossover from Bose-Einstein condensation (BEC) 
to Bardeen-Cooper-SchriefFer (BCS) state has drawn re- 
newed interest in past few years due to rapid experimen- 
tal progress in cold atomic two-cornponent Fermi gases. 
Several experimental groups [l|, 0, H, 3, H, 0, Q achieved 
the BEC-BCS crossover in the two-component atomic 
Fermi gases due to the magnetized Feshbach resonance 
mechanism [HI. The fermionic system becomes 

molecular BEC for strong repulsive interaction and trans- 
form into the BCS states when the interaction is attrac- 
tive. Near the resonance, the zero energy s-wave scat- 
tering length a exceeds the interparticle spacing and the 
interparticle interactions are unitarity limited and uni- 
versal. 

As in the case of bosonic clouds the frequencies of col- 
lective modes of Fermi gases can be measured to high ac- 
curacy. The collective oscillation frequencies of a trapped 
gas can provide crucial information on the equation of 
state of the system. The experimental results on the col- 
lective frequencies of the lowest axial and radial breath- 
ing modes on ultra cold gases of ®Li across the Feshbach 
resonance have also become available [H, • Since the 
weak-coupling BCS and unitarity limits are character- 
ized by the same collective oscillation frequency, it is an 
interesting to find out another observable which makes 
a clear identification of these two regimes and to better 
characterize two kinds of superfluid. 

All the experiments of cold atomic Fermi gases are 
done in a cigar shaped geometry in which the atomic 
density is inhomogeneous in the radial plane and quasi- 
homogeneous along the symmetry axis. The axial ex- 
citations of a cigar shaped Fermi superfluid can be di- 
vided into two regimes: short wavelength excitations 
whose wavelength is much smaller than the axial size and 
long wavelength excitations whose wavelength is equal 
or larger than the axial size of the system. In the later 
case, the axial excitations are discrete and the lowest ax- 
ial breathing mode frequency has been measured p^ . 
In the former case, the short wavelength axial phonons 
with different number of radial modes of a cigar-shaped 



Fermi superfluid give rise to the multibranch spectrum 
. These are similar to the electromagnetic wave prop- 
agation in a wave guide. In the BCS side of the resonance, 
the low energy bosonic excitations, apart from the gaped 
fermionic excitations, of the Fermi superfluid are called 
multibranch Bogoliubov- Anderson (BA) modes. In the 
usual electronic superconductors, the BA phonon mode 
is absent due to the long-range Coulomb interaction. 

In this work we find that the low energy multibranch 
modes in the BCS limit are different from that of in 
the unitarity limit, although discrete radial and axial 
mode frequencies are same due to the same exponent in 
the equation of state in the unitarity and BCS regimes. 
Therefore, these multibranch modes can be used to iden- 
tify and characterize the weak-coupling BCS and unitar- 
ity states. 

It is an interesting to study how one can probe such 
bosonic modes in the current available experimental 
setup. In fact, these multibranch low-energy bosonic 
modes could be observed by measuring dynamic struc- 
ture factor (DSF) in a Bragg scattering experiment. 
Bragg spectroscopy of a trapped atomic system has 
proven to be an important tool for probing many bulk 
properties such as dynamic structure factor p^ . verifi- 
cation of multibranch Bogoliubov excitation spectrum in 
the usual atomic BEC [16], correlation functions and mo- 
mentum distributions of a phase fluctuating Bose gases 
[itL [TsI I . There has been a suggestion of Bragg scatter- 
ing experiment to study the fermionic excitations at zero 
temperature [l9j as well as at temperature close to the 
critical temperature [lO]. 

In this work, we also calculate DSF in the various 
regimes of the crossover, namely the weak-coupling BCS, 
unitarity and molecular BEC regimes. In actual Bragg 
scattering experiments, the measured response function 
is momentum transferred Pz{t) to the superfluid by the 
Bragg pulses. The momentum transferred is directly re- 
lated to the dynamic structure factor when time duration 
of the Bragg pulses is long enough. Therefore, we also 
study momentum transferred to the superfluid by the 
Bragg pulses. In order to probe the multibranch nature 
of the modes, we also estimate required values of wave 
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vector and time duration of the two-photon Bragg pulses 
in the Bragg scattering experiments. 

Recently, there is a measurement of sound velocity 
along the crossover '2l'| . The measured sound velocity do 
not match with the theoretical prediction in the molecu- 
lar BEC regime, however it matches very well in other 
regimes of the crossover [22|. One can also estimate 



the sound velocity in the three different regimes along 
the crossover by measuring slope of the phonon mode. 
It may resolve the puzzle of mismatch of the measured 
sound velocity with that of the theoretical prediction in 
the molecular BEC state. 

This paper is organized as follows. In section II, we 
provide quantized hydrodynamic description of the Fermi 
superfluid along the crossover. In section III, we calculate 
the dynamic structure factor along the crossover. In sec- 
tion IV, we discuss the possible Bragg scattering experi- 
ment in this system and present results of the momentum 
transferred to the system due to the Bragg pulses. We 
also discuss a summary and conclusions in section V. 



II. QUANTIZED HYDRODYNAMIC THEORY 
OF FERMI SUPERFLUID 

We use hydrodynamic model with the Weizsacker 
quantum pressure term to describe low-energy dynam- 
ics of a two-component Fermi superfluid at zero temper- 
ature. This system can be well described by a time- 
dependent non-linear Schrodinger equation as follows 



ih 



dt 



2M 



V2-|-K,,t(r,z)-f /2(n)]V-(r,i), (1) 



where the non-linear term JjL{n) = ^[ne(n)] is the chem- 
ical potential in a uniform system and e(n) is the ground 
state energy per particle. Here, M is the mass of a Fermi 
atom and Vcxt('', z) = (Af/2)(a;^r^ -l-cj^z^) is an external 
harmonic trap potential with lo^ >> oJz- We have taken 
u)z/^^>r = 0.1 and total number of atoms N = 2.0 x 10^ 
in all numerical calculations. On the basis of quantum 
Monte Carlo data of Astrakharchik et al. Manini 
and Salasnich [2^ proposed a very useful analytical fit- 
ting expression for e(n) (see Eq. (6) and Table 1 of Ref. 
[25j). The Hamiltonian corresponds to Eq. ([1]) can be 
written as 

H = J drr{r,t)[-—^^ + V,^t{r,z)+fl{nmr,t). 

Using phase (6')-density (n) representation of the or- 
der parameter of the composite bosons: '0(r, = 
\/n{r, t)e^^^'^'*\ the above Hamiltonian becomes, 



h 1 
H = I dr[\^{-—V^\^) + -Mv^n + Vcy,t{r,z)h 



Linearizing the density and phase around their equi- 
librium values: n(r, i) = no{r) + 6n{r,t) and v(r, i) — 
(5v(r,t), where v = {h/My^O is the superfluid velocity. 
Keeping upto quadratic fluctuations, the above Hamilto- 
nian reads 



H = Hn 



dr 



(3) 



where Hq is the ground state energy. By using time- 
dependent Heisenberg equations of motion for the den- 
sity and velocity fluctuations, one can get continuity and 
Euler's equations which are given by 



and 



Sn = —V • [no(r)(5v] 



M6^ ^ -V[^lJn]. 



(4) 



(5) 



Since the Hamiltonian is quadratic in terms of the fluc- 
tuation operators, it can be diagonalized by using the 
following standard canonical transformations: 

Shir,z,t) - ^[A,-fc^,-fc(r)d,-fee^('=^-"^W*) + /i.e.], (6) 



and 



(r, z, t) = Y,[Bj.k^j.k{r)aj.ke'^'''-^^^''^'^ + h.c]. (7) 



Here, j is a set of two quantum numbers: radial quan- 
tum number, rir and the angular quantum number, m. 
Also, k is the axial wave vector. The density and phase 
fluctuations satisfy the following equal-time commutator 
relation: [5n{r), 50{r')\ — iS{r — r'). One can easily show 
that 



A 



3,k 



I TiLOj (k) 



B 



j-k 



flLOj{k) 



(8) 



(2) 



We are assuming that tpj^kir) satisfies the orthonor- 
mal conditions: / drip* f.{r)ipji ^(r) = djj> and 

E,^.,fcWV'*fe(r') = '5(r-r'). 

We assume power-law form of the equation of state as 
/i(n) — CrP . Here, C depends on interaction strength 
and the effective polytropic index 7 is a function of a 
dimensionless parameter y = l/kpa, where kp is the 
Fermi wave vector and a is the scattering length between 
Fermi atoms of different components. The weak-coupling 
BCS [y « —1) and the unitarity {y = 0) states are 
described by the same exponent 7 = 2/3 with differ- 
ent values of C. For the weak-coupling BCS regime, 
Cbcs ^ {3Tr'^)'^/^{3n'^/10M) and for the unitarity regime, 
Cuni = 0.44Cbcs- The molecular BEC state (y >> 1) 
is described by 7 = 1 and Cbcc = Anh'^am/M , where 
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flm = 0.6a is the molecular scattering length [26[. In 
our calculation we have assumed a™ = 1.0 x 10~®to. 
The power-law form of the equation of state is being 
used successfully to study the Fermi superfluid along the 
crossover [25i, i27|, i28j, i2^, ,301 ■ equilibrium, the density 
profile takes the form no(r) = (/i/C)^/'''(l — f^)^/'', where 
f = r/Ro, Ro = y/2n/Muj^ and fj, is the chemical poten- 
tial in the non-uniform system, which can be obtained 
from the normalization condition. 

Taking first-order time-derivative of Eq. (jj]) and using 
Eq. ([5]), the second-order equation of motion for the 
density fluctuation is given by 



= V - [no(r)V-g^|„=„„(5n]. 



(9) 



Using the polytropic form of the equation of state and 
Eq. ([6]), then Eq. ([9]) reduces to the following equation: 



(10) 



where a) = uj/tOr and k — kRQ. 

For /c = 0, it reduces to a two-dimensional eigenvalue 
problem and the solutions of it can be obtained analyti- 
cally. The energy spectrum is given by 



w- = \m\ + 2nr[^{nr + |m|) + 1]. 



(11) 



The corresponding orthogonal eigenfunction is given by 

V^j- CX (1 - f2)l/7-l^;|™lpa/7-l,l™l)(2f2 „ ly^'t'^ (12) 

where P^'^\x) is a Jacobi polynomial of order n and (j) 
is the polar angle. 

The solution of Eq. (fTU]) can be obtained for arbitrary 
value of k by numerical diagonalization. For fc 7^ 0, we 
expand the density fluctuation as 



(13) 



Substituting the above expansion into Eq. (jlOp . we 
obtain. 



be obtained separately in each angular momentum sub- 
space. We can obtain all low energy multibranch spec- 
trum on the both sides of the Feshbach resonance includ- 
ing the unitarity limit from Eq. (HH). Equations 
and (jisp show that the spectrum depends on the average 
over the radial coordinate and the coupling between the 
axial mode and transverse modes within a given angular 
momentum symmetry. Particularly, the coupling is im- 
portant for large values of k. We are interested to study 
m = states since these states are excited in the Bragg 
scattering experiments due to axial symmetry of the sys- 
tem. We show low energy multibranch modes of m = 
states in three different regimes in Fig. 1. The top three 
panels of Fig. 1 show the multibranch spectrum in the 
three different regimes. To compare these spectrum, we 
have plotted all those spectrum of different regimes in a 
single frame, which is shown in the bottom panel of Fig. 
1. The lowest branch corresponds to the Bogoliubov ax- 
ial mode with no radial nodes. This mode has the usual 
form Lj — Cgk aX low momenta, where Cg is the sound 
velocity. In the limit of small k, the other branches have 
free-particle dispersion due to the gaped nature of these 
modes. 

The discrete radial and axial modes are same in the 
unitarity and BCS regimes since the exponent 7 are the 
same for both the regimes. However, the multibranch 
modes are different in the unitarity and BCS regimes in 
spite of the same exponent in the equation of state. This 
is due to the different radial sizes in those two regimes for 
a given number of atoms and the trap potential. There- 
fore, these low energy axial propagation of discrete radial 
modes can be used to characterize different regimes of the 
superfluid. The density fluctuations for a fixed value of 
the axial momentum are plotted in Fig. 2. The den- 
sity fluctuations corresponds to the multibranch modes 
are also different in the different regimes. In Fig. 2, 
the magnitude of the density fluctuations are given in an 
arbitrary unit since these are the linear fluctuations. 



III. DYNAMIC STRUCTURE FACTOR 

The dynamic structure factor is the Fourier transfor- 
mation of density-density correlation functions and it is 
given as 



= [Cjj{kf -[\m\+2nr{-i{nr + \m\) + l)] 



(14) 



j' 



Here, the matrix element Mjji is given by 

A/„-, = / d''rii)if^il}i'. (15) 



The above eigenvalue problem (Eq. (fT4|) ) is block diag- 
onal with no overlap between the subspaces of different 
angular momentum, so that the solutions to Eq. (|14p can 



S{k,uj) = / dvdv'dt < Sn\r,t)Sh{r',0) > e'^' 



This can be written as 

5(fc,c.) = ^5,(fc)<5(c.-^,(fc)), (16) 



where the weight factor Sj{k) is given by 



S,{k) 



hujj{k) 



(17) 
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FIG. 1: (Color online) Plots of the low energy multibranch 
modes in the three diflerent regimes. 



Here, — J dre~^^''^ipj^k{r) is the Fourier transform 

of the eigenfunctions 4'j,k{f')- The weight factors are plot- 
ted in Fig. 3. The weight factors Sj{k) determine how 
many modes are excited for a given value of k. For ex- 
ample, when kar = 1.0, = {0, 1}, rir = {0, 1, 2 } 
and rir = {0, 1, 2, 3} modes are excited in the molecu- 



lar BEC, unitarity and weak-coupling BCS regimes, re- 
spectively. The harmonic oscillator length is defined as 
Or = ^jTi/MuJr. To excite many other low-energy modes, 
the wave vector in the Bragg potential must be large. 
From Fig. 3, it is clear that S'o(fc) ^ k and 5'„^>o(fc) 
is almost constant in the limit of small k. Therefore, 
^o(k) ~ k and LOn^-yQ{k) ^ k^ when k is very small. 
Therefore, our analysis also satisfies the Feynman-like 
relation u;„^(/c) ~ k'^/Sn^{k). The dynamic structure 
factors of the three different regimes for ka^. = 0.5 are 
plotted in Fig. 4. The delta function in Eq. (|16p is 
replaced by the Lorcntzian form to plot Fig. 4. 
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FIG. 2: (Color online) Plots of the low-energy density fluctu- 
ations in the three different regimes for kar = 0.5. 



5 



0.8 



0.6 



0.4 



0.2 




0.4 



0.2 






1 


2 

q 


3 4 


5 
















































Unitarity 




D 


1 


2 

q 


3 4 


5 



0.8 r 






0.7 


/' \ ,- M. 




0.6 


/ II. 3 




0.5 - 






0.4 - 


/ III nr=^^^^^^^^n>. 




0.3 






0.2 




Weak-coupling BCS 


0.1 1 






L 







q 



FIG. 3: (Color online) Plots of the weight factors in the three 
different regimes. 



IV. BRAGG SCATTERING EXPERIMENT 

The behavior of these multiple peaks in the dynamic 
structure factor can be resolved in a two-photon Bragg 
spectroscopy, as shown by Steinhauer ei al. [Slj for 
usual BEC. In the two-photon Bragg spectroscopy, the 
dynamic structure factor can not be measured directly. 
Actually, the observable in the Bragg scattering exper- 
iments is the momentum transferred to the superfluid, 
which is related to the dynamic structure factor and re- 
flects the behavior of the quasiparticle energy spectrum. 
The populations in the quasiparticle states can be con- 
trolled by using the two-photon Bragg pulse. When the 
superfluid is irradiated by an external moving optical po- 
tential V^p — Vsit) cos{qz — ojt), the excited states are 
populated by the quasiparticle with energy hco and the 
momentum hq, depending on the value of q and lo of 
the optical potential Vop- Here, Vb is the intensity of 
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FIG. 4; (Color online) Plots of the dynamic structure factors 
of the three different regimes for kar = 0.5. 



the Bragg pulse. Suppose the system is subjected to a 
time-dependent Bragg pulse which is switched on at time 
t > and q is also along the z-direction. We calculate, 
similar to the calculation of Refs. [13, [s^l , the momen- 
tum transfer to the superfluid from the moving optical 
potential and it is given by 



2h 



X J2hqS,iq)x[F,{Lu^t)-F,{LU+t)], (18) 
j 

where aj^k{t) is the time-evolution of the quasiparticle 
operator of energy fiiOjik) and 

n[(uj,{q)±u)tl2] V 
(c.,(g)±u;)/2 ) 



(19) 



For positive lo and a given q such that Sj (q) is maximum, 
the momentum transferred Pz{t) is resonant at the fre- 
quencies to — t-Ojiq)- The width of the each peak goes 
like 27r/t. For large t and uj^ << Ur, one can show that 

In Fig. 5, we plot the net momentum transfer Pz{t) vs 
the Bragg frequency ui for three different choices of the 
time duration of the Bragg pulses. Figure 5 shows that 
the shape of the Pz{t) strongly depends on the time du- 
ration of the Bragg pulses Ib- When tB — Q-5Tr, Pz{t) 
is a smooth curve with a single peak. Here, we define 
Tr — 27r/ti)j, is radial trapping period. When Ib = l.OT^, 
there is a little evidence of few small peaks start devel- 
oping in Pz{t)- When Ib = 2.0Tr, the multiple peaks 
in Pz{t) appears prominently. Therefore, the duration 
of the Bragg pulses t b should be greater than the radial 
trapping period in order to resolve different peaks in 
the DSF. Figure 5 also shows that the number of quasi- 
particle modes in unitarity regime is much higher than 
the other regimes of the crossover. This is due to the fact 
that Pzit) is proportional to the number of quasiparticle 
modes for a given k. 
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FIG. 5: (Color online) Plots of the momentum transferred 
in the three different regimes for fca^ = 0.5 and for different 
time duration of the Bragg pulses: t = O.STr (solid), t — l.OTr 
(dashed) and t = 2.0Tr (dot-dashed). 



SUMMARY AND CONCLUSIONS 



We have presented the quantized hydrodynamic the- 
ory of cigar shaped Fermi superfLuid along the BEC- 
BCS crossover by using the power law form of the equa- 
tion of state. We have calculated multibranch low en- 
ergy bosonic modes and the corresponding density fluc- 
tuations in three different regimes along the BEC-BCS 
crossover, namely the weak-coupling BCS, unitarity and 
molecular BEC states. Then we have presented results of 
the dynamic structure factor calculation. We have also 
calculated the momentum transferred to the superfluid 
by the Bragg pulses and shown that the multibranch na- 
ture would be observed when the time duration of the 
Bragg pulses is greater than the radial trapping period. 

We have found that the axial propagation of discrete 
radial modes in the weak-coupling BCS and unitarity 
regimes are different, although the axial and radial modes 
are same in both cases. One can identify these two 
regimes by probing these multibranch modes with the 
help of the Bragg scattering experiment. We have also 
seen that the number of quasi particle modes in the uni- 
tarity regime is quite large compared to other two regimes 
for a fixed value of the Bragg momentum and the number 
of atoms. Therefore, the response function in the unitar- 
ity regime will be more prominent than the other two 
regimes in the crossover. Moreover, one can estimate the 
sound velocities in different regimes along the crossover 
by measuring the slope of the phonon modes. 
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